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Abstract 

We study anode-directed ionization fronts in curved geometries. When the magnetic effects 
can be neglected, an electric shielding factor determines the behavior of the electric field and the 
charged particle densities. From a minimal streamer model, a Burgers type equation which governs 
the dynamics of the electric shielding factor is obtained. A Lagrangian formulation is then derived 
to analyze the ionization fronts. Power laws for the velocity and the amplitude of streamer fronts 
are observed numerically and calculated analytically by using the shielding factor formulation. 
The phenomenon of geometrical diffusion is explained and clarified, and a universal self-similar 
asymptotic behavior is derived. 

PACS numbers: 52.80.Hc, 05.45.-a, 47.54.+r, 51.50.+V 
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I. INTRODUCTION 



In a perfect dielectric medium, charged particles form electrically neutral atoms and 
molecules due to powerful electric forces. Since there are not free charges in this medium, 
electric current does not flow inside it. However, if a very strong electric field is applied to 
a medium of low conductivity in such a way that some electrons or ions are created, then 
mobile charges can generate an avalanche of more charges by impact ionization, so that a 
low temperature plasma is created, and an electric discharge develops Q. This process is 
called electric breakdown and it is a threshold process: there are not changes in the state 
of the medium while the electric field across a discharge gap is gradually increased, but at 
a certain value of the field a current is created and observed. 

A streamer is a ionization wave propagating inside a non-ionized medium, that leaves 



a non-equilibrium plasma behind it. They appear in nature and in technology 



a a. 



A 



streamer discharge can be modelled using a fluid approximation based on kinetic theory (^J . 
Defining the electron density N e as the integral of the electron distribution function over all 
possible velocities, we get 

dN 

^ + V R -J e = S e , (1) 

where r is the physical time, Vr is the gradient in configuration space, U e (R, r) is the 
average (fluid) velocity of electrons, S e is the source term, i.e. the net creation rate of 
electrons per unit volume as a result of collisions, and J e (R, r) = iV e (R, r) U e (R, r) is the 
electron current density. Similar expressions can be obtained for positive N p and negative 
N n ion densities. 

A usual procedure is to approximate the electron current J e as the sum of a drift (electric 
force) and a diffusion term 

J e = -fi e £N e - D e V R N e , (2) 

where £ is the total electric field (the sum of the external electric field applied to initiate the 
propagation of a ionization wave and the electric field created by the local point charges) 
and n e and D e are the mobility and diffusion coefficients of the electrons. Note that, as the 
initial charge density is low and there is no applied magnetic field, the magnetic effects in 
equation (J2J) are neglected. This could not be done in cases where the medium is almost 
completely ionized or cases in which an external magnetic field is applied, leading to different 
treatments 



Some physical processes can be considered giving rise to the source terms. The most 
important of them are impact ionization (an accelerated electron collides with a neutral 
molecule and ionizes it), attachment (an electron may become attached when collides with a 
neutral gas atom or molecule, forming a negative ion), recombination (a free electron with a 
positive ion or a negative ion with a positive ion) and photoionization (the photons created 
by recombination or scattering processes can interact with a neutral atom or molecule, 
producing a free electron and a positive ion) . 

It is also necessary to impose equations for the evolution of the electric field S. It is usual 
to consider that this evolution is given by Poisson's law, 

V R • S = - (N p - N n - N e ) , (3) 

where e is the absolute value of the electron charge, Eq is the permittivity of the gas, and 
we are assuming that the absolute value of the charge of positive and negative ions is e. 

Some simplifications can be made when the streamer development out of a macroscopic 
initial ionization seed is considered in a non-attaching gas like argon or nitrogen. For such 
gases, attachment, recombination and photoionization processes are usually neglected. A 



minima 



model turns out and has been used to study the basics of streamer dynamics 
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ll| . In those cases the evolution of electron and positive ion densities in early 
stages of the discharge can be written as 

-gf = Vr ■ (/i,£ N e + D e V R N,) + v,N c , (4) 

On time scales of interest the ion current is more than two orders of magnitude smaller than 
the electron one so it is neglected in (jSJ). In these equations ViN e is a term accounting for 
impact ionization, in which the ionization coefficient z/j is given by the phenomenological 
Townsend's approximation, 

Vi = fi e \£\a e- £o/l£ \ (6) 

where /i e is the electron mobility, «o is the inverse of ionization length, and Sq is the char- 
acteristic impact ionization electric field. 

It is convenient to reduce the equations to dimensionless form. The natural units are 
given by the ionization length R = a^ 1 , the characteristic impact ionization field £q, and 
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the electron mobility fi e , which lead to the velocity scale Uo = fi e £o, and the time scale 
t = Rq/Uq. The values for these quantities for nitrogen at normal conditions are a^ 1 ~ 
2.3 fim, £ ~ 200 kV/m, and \i e ~ 380cm 2 /(Vs). We introduce the dimensionless variables 
r = H/Rq, t = t/tq, the dimensionless field E = £/Sq, the dimensionless electron and 
positive ion particle densities n e = N e /N and n p = N p /N with N = e £ /(eRo), and the 
dimensionless diffusion constant D = D e / '(RqIIq). 

In terms of the dimensionless variables, the minimal model equations become 

^ = V-j + nJ(|E|), (7) 

^ = nJ(|E|), (8) 

n p - n e = V ■ E, (9) 

j = n e E + D Vn e , (10) 

/(|E|) = lEle- 1 ^, (11) 
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for the min- 



where V = V r , and j is the dimensionless electron current density. 

Some properties of planar fronts have been obtained analytically |6 
imal model. A spontaneous branching of the streamers from numerical simulations has 
been observed 0, as it occurs in experimental situations Q|. In order to understand this 
branching, the dispersion relation for transversal Fourier-modes of planar negative shock 
fronts (without diffusion) has been derived ^(|. For perturbations of small wave number k, 
the planar shock front becomes unstable with a linear growth rate l-E^I^- It has been also 
shown that all the modes with large enough wave number k (small wave length perturba- 
tions) grow at the same rate (it does not depend on k when k is large). However, it could 
be expected from the physics of the problem that a particular mode would be selected. To 
address this problem, a possibility is to consider the effect of diffusion. It is also interesting 
to investigate the effect of electric screening since, in the case of curved geometries, this 
screening might be sufficient to select one particular mode. 

In this paper, we study the properties and structure of anode-directed ionization fronts 
with zero diffusion coefficient for curved geometries. We start discussing the consequences 
of neglecting the magnetic field effects in the physics of the streamer evolution. As a con- 
sequence, an electric shielding factor can be introduced which determines the behavior of 
the electric field and the particle densities. From the minimal streamer model, a Burgers 
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type equation which govern the dynamics of the electric shielding factor is deduced. This 
allows us to consider a Lagrangian formulation of the problem simplifying the analytical and 
numerical study of the fronts. We apply this new formulation to planar as well as curved 
geometries (typical in experimental set-ups). Power laws for the velocity and the amplitude 
of streamer fronts are observed numerically. Theses laws are also calculated analytically by 
using the shielding factor formulation. The geometrical diffusion phenomenon presented in 
[l5l | is explained and clarified, and a universal self-similar asymptotic behavior is derived. 

The organization of the paper is as follows. In Section II, we will show that all the physics 
involved in the minimal model can be rewritten in terms of the electric shielding factor, that 
determines the behavior of the charge densities and the local electric field in the medium. 
This allows a simple analysis of the model when written in Lagrangian coordinates. Within 
this framework, we perform in Section III, as an illustration of the Lagrangian formulation, 
the analysis of planar fronts (without diffusion). In Section IV, we study the evolution of 
ionization fronts in which the initial seed of ionization is such that the electron density 
vanishes strictly beyond a certain point for cylindrical and spherical symmetries. We obtain 
precise power laws for both the velocity of the moving fronts and their amplitude. In Section 
V, we analyze the special features that appear if the initial seed of ionization is not completely 
localized but the charge densities slowly decrease along the direction of propagation. For 
curved geometries, this initial distribution gives rise to a new diffusion-type behavior that we 
call geometrical diffusion. A universal self-similar asymptotic shape of the fronts is predicted 
and observed. In Section VI, we establish our conclusions. 



II. ELECTRIC SHIELDING FACTOR 



In this section we will reformulate the problem of the evolution of streamer fronts in 
he minimal model by introducing a new quantity called the electric shielding factor, as in 
I5I 1 . The equation describing the evolution of the shielding factor makes easier the study of 
curved ionization fronts. 

We begin with a brief discussion about the consequences of neglecting the magnetic effects 
in the minimal streamer model. In this model, it is assumed that the magnetic field effects 
are negligible, in a first approximation, because (i) the fluid velocity of the electrons is 
much smaller than the velocity of light, and (ii) the initial magnetic field is zero. Strictly 
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speaking, if the magnetic field is zero in the evolution of the ionization wave, then Faraday's 
law implies that the electric field is conservative (i.e. V x E = 0). This means that, in cases 
in which the evolution of the ionization wave is symmetric (planar, cylindrical or spherical), 
the electric field would evolve according to this symmetry, so one can write 



where Eo(r) is the initial electric field (that is conservative since it is created by an applied 
potential difference) and u(r, t) is some scalar function with the same symmetry as the initial 
electric field, since 



and therefore V« is parallel to E . Consequently, the relation (J 12)) assures that the magnetic 
field will always be zero if the initial magnetic field is zero and it is a direct consequence of 
the hypothesis of the minimal streamer model. However, experimental observations indicate 
that streamers can change their direction while evolving, suggesting that the ionization 
process can be non-symmetric in some cases and that the local magnetic field may play a 
role, pointing to future modifications of the model. Nevertheless, as a first approach to the 
problem, we will consider here symmetric situations. In these cases, the minimal model can 
be applied for the streamer evolution, and the relation (JT2jl is strictly correct. 

The above discussion leads to an unexpected consequence on the minimal streamer model: 
the quantity u defined in f|12|) determines completely, without any physical approximation, 
the electric field and the particle densities during the evolution of the ionization wave if 
the diffusion is neglected. If we take the diffusion coefficient D equal to zero, the minimal 
streamer model given by equations ([7f)-([TT|) can be rewritten as 




(12) 



= V x E = E x Vw, 



(13) 




(14) 



drip 



n e |E|e 



1/|E| 



(15) 
(16) 



dt 
V- E 



Subtracting equation (JT3J) from ((T3J) 



we obtain 




(17) 



By taking the time derivative in equation (|16|) . we obtain 




(18) 
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V- (f^ + nJZ) =0. (19) 



and hence, using (JUj), we get 



\dt 

Since the electric current is given by n e E, expression (fT^j) is simply the divergence of 
Ampere's law applied to our case, with the right hand side being the divergence of the 
curl of the magnetic field, which is always zero. In the particular situation in which the curl 
of the magnetic field in the gas is negligible, as it occurs in the framework of the minimal 
model for symmetric situations (as discussed above), this expression can be also written as 

— + n e E = 0. (20) 

This is a linear first-order ordinary differential equation for the electric field, so that it can 
be trivially integrated to give 

E(r, t) = E (r) exp (- J rfrn e (r, r)\ , (21) 

which supplies the local electric field E in terms of the initial electric field E multiplied 
by the electron density n e integrated in time. The physical behavior of the electric field 
screened by a charge distribution suggests that an important new quantity can be defined 

as 

u(r,t) = exp ^- J drn e (r,r)Y (22) 

so that the relation ()12j) is re-obtained. This means that, when u is determined in a particular 
situation, the electric field E is known. Moreover, using equations ()22|) and (jlfijl . we obtain 
that the particle densities are also determined by u and the initial condition Eo(r) for the 
electric field, through 

i A 1 du{r,t) 

ne(M) = - U (r,t) dt ' (23) 

n P (r, t) = - ^^^f + V ■ (Eo(r)n(r, *)) . (24) 

Equation ()12j) reveals clearly the physical role played by the function u(r,t) as a factor 
modulating the electric field E(r, t) at any time. For this reason, u can be termed shielding 
factor and determines a screening length that depends on time. This is a kind of Debye's 
length which moves with the front and leaves neutral plasma behind it Q]. As the shielding 
factor determines the particle densities, equation ()13|) implies that the particle densities have 
the same symmetry as the initial electric field. 



The definition of the shielding factor u and the mathematical treatment explained above 
reduce the problem of evolution of charged particle densities and electric field in the gas to a 
simpler one: to find equations and conditions for the shielding factor u(r,t) from equations 
and conditions for the quantities E, n e and n p . Substituting equations (|12j) - (j24j) into the 
original model equation f|14| ) — f)lfc>| ) . we find 



d ( 1 du 



dt \u dt 



- V ■ (E u 



' ° l dt 



(25) 



where |Eo| is the modulus of the initial electric field Eo- The last term in this expression 
can be written as 



, du 



|E |— eTO^ 



— / e- 1/s ds 
dt J 



so that 



_ ___ V .(E 0M )U-/ ,-,/, 



This equation can be integrated once in time to give 



1 du 



-rZ-V.(E «)= / e-V'ds + Gir), 
u dt J Q 



(26) 



(27) 



(28) 



where the function G(r) is given by 



E \u 



- 1/s ds 



(29) 



t=0 



The initial conditions for u and du/dt can be easily related to initial conditions for particle 
densities using (IT2*|) and The results are 



9m 

at 



(30) 
(31) 



t=0 



where n e0 (r) is the initial value of n e (r, t). Then, 



G*(r) = -n e0 (r)-V-E - 



|Eo| 



(32) 



which, if n p o(r) is the initial value of the dimensionless ion density, can also be written as 

|Eo| 



G(t) = -n pQ (r) - / e-V'ds. 
Jo 



(33) 
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As a consequence, the evolution of u(r,t) is given by 




l/s 



ds, 



(34) 



u(r,0) = u (t) = 1, 



(35) 



with appropriate boundary conditions which should be imposed depending on the particular 
physical situations one wishes to consider. Note that we have written the complete minimal 
model (with D = 0) in one single equation for the shielding factor u. All the physics in the 
minimal model is contained in the evolution equation (|34|). The shielding factor is related 
to charged particle densities and electric field through expressions (H2J), (|23Jl and (J2U). This 
formulation ()34|) allows us a much simpler analysis than the original one, and will provide 
us with some insight into some unsolved problems on streamer formation. 

III. THE PLANAR CASE AND THE LAGRANGIAN FORMULATION 

In this section, we will use the shielding factor in order to find the main features of planar 
anode-directed ionization fronts without diffusion. This is a very simple way of testing the 
usefulness of the new formulation for further generalization, as the results can be compared 
with that found in [6] using a different approach. 

A. Lagrangian formulation 

We consider an initial experimental situation as follows. Two infinite planar plates are 
situated at z = and z = d respectively (z is the vertical axis). The space between the plates 
is filled with a non-attaching gas like Nitrogen. A stationary electric potential difference Vq 
is applied to these plates, so that V(d) — V(0) = Vq > 0. To initiate the avalanche, an initial 
neutral seed of ionization is set at the cathode, so that n e0 (z) = n p o(z) = po(z)- We study 
the evolution of negative ionization fronts towards the anode at z — d. 

As the applied potential is constant, the initial electric field Eo between the plates results 

in 



Eo = — E u z , E = —. 

d 

It is useful for the computations to define the coordinate x as 



(36) 



z 



(37) 
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so that the evolution of the shielding factor u (|34|) is given as the solution of the equations 



This is a typical Burgers type equation with an integral term. Hence, we can use some of 
the classical techniques developed to deal with this equation. In particular, we can integrate 
along characteristics and transform equation (JHS|) into the system 



which yields a Lagrangian formulation of the problem. The solutions of this dynamical 
system with initial data given by x(0) = x , u(0) = 1 allows us to compute the profiles 
for u(x,t) at any time. Then using equations (fT2|) . (f^|) and (j2U), it is possible to trace 
the profiles of the electric field or the charge densities at different times. This has been 
done in Fig. ^ in which the electron density n e is plotted as a function of the coordinate 
x = z/Eq. We have chosen a neutral initial seed of ionization sufficiently localized near the 
negative plate, i.e. the electron and positive ion densities are initially equal and, moreover, 
they vanish beyond a certain point in the x axis (mathematically, this situation is described 
by saying that the initial condition is of compact support). After evolution, the electron 
density converges to a travelling wave, as can be seen in Fig. [0 This travelling wave has a 
constant propagation velocity and a constant amplitude, as can be seen in the figure, and it 
is a shock front. This shock front appears only if the initial condition is of compact support, 
as we will see later. 

B. Analytical computations 

The fact that, in the planar case, the integral term in equation ()38|) does not depend 
explicitly on x has an interesting consequence: the velocity of the front can be computed 
directly from the equation for the shielding factor and it is completely determined by the 
initial condition. Let us assume that the initial density po(x) (for both electron and positive 
ion densities) decays sufficiently fast as x goes to infinity. Then we can neglect the term 




(38) 
(39) 



dx 



(40) 




(41) 
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FIG. 1: Electron density n e of a planar ionization wave for fixed time intervals vs coordinate 
x = z/Eq, in which z is the propagation direction (from the negative to the positive plate) and Eq 
is the modulus of the initial electric field (that is constant) between the plates. The initial data 
is a compactly supported neutral seed of ionization near the negative plate x = 0. By using the 
formulation in terms of the shielding factor we can see that, after evolution, the ionization wave 
converges to a travelling wave of constant velocity and constant amplitude, called a shock front. 

upo(x) in (J38|) for x ^> 1 and look for a solution of the resulting equation in the form 

u(x,t) = f{£ = x-ct). (42) 

Using the equation (|38|) for the evolution of the shielding factor with this approximations, 
we obtain the differential equation 

df Mf) 



di c-f 

where the quantity v(f) is given by 



(43) 



fE 

v (f) = / e - 1/s ds. (44) 

The front of the wave is localized, at a given time, in points in which / — 1. In these points, 
up to first order in /, the quantity v(f) results in v(f) ~ EQe~ l l E °(l — f). Inserting these 
approximations into equation we get 

df Er,e~ 1/E ° 



By integrating this expression and using (|23|) to obtain the electron density, it can be easily 
seen that physically acceptable solutions of this equation (i.e. positive value of the electron 

11 



density in all points), correspond only to values of / given by / < 1 for large £. From (|45jh 
these physical solutions appear only if c > 1 . So that the velocity c z of propagation of the 
front, in the original z coordinate, satisfies 

c z > E , (46) 

in agreement with the result found in 

Moreover, by using this formulation it is also possible to link the asymptotic behavior 
of the initial condition n e0 with the propagation velocity, so that it will be shown that the 
initial condition determines the velocity of the front. Suppose that the initial condition for 
the electron density behaves like n e o ~ Ae~ Xx as x —* oo. Then, the asymptotic behavior of 
the travelling wave satisfies 

n e ~ Ae~ x ^ , as £ -> oo. (47) 
Using the relation (|2Hj) . this means that the shielding factor u = /(£) behaves like 

/ ~ 1 - ^e" A « , as £ ^ oo. (48) 

When this expression is introduced into equation (|45p. the relation 

A = EoeX c ^- 1/Eo) , (49) 

appears. This is the way in which the asymptotic behavior of the initial condition determines 
the propagation velocity. By using this link, in Fig. |2]we have plotted several travelling wave 
profiles for different values of c. In Fig. EJleft), the shielding factor has been plotted, and in 
Fig. fright), the corresponding electron density, both as a function of £ = x — ct. 

As mentioned above, the minimum value of the propagation velocity c of the travelling 
wave is c = 1 (it corresponds to c z = Eq). For this velocity, a shock front can be clearly 
seen in Fig. El Using equation (JH2J), this shock front appears when A — >• oo, so that the 
initial condition is of compact support, i.e. the initial distribution of charge vanishes strictly 
beyond a certain point. In Fig. Q an initial condition fulfilling these requirements has been 
chosen and a shock front has appeared as predicted by (|4*9*|) . The amplitude of the shock 
front can be also obtained from equation (|45jl and is given by 

n e ~E e- 1/Eo . (50) 

The rest of profiles in Fig. El correspond to travelling waves in which the velocity c is larger 
than 1, so the initial distribution does not have compact support. 
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FIG. 2: (Left) Electric shielding factor u, and (Right) electron density n e of several travelling wave 
profiles. In the figures, £ = x — ct and c takes values from c = 1 (it corresponds to a shock front) to 
c = 6. The propagation velocity depends on the way in which the initial electron density behaves 
as x goes to oo. 

C. Accelerated fronts 



The treatment given above could suggest that ionization fronts always move with constant 
velocity if the physical setting has planar symmetry. However, it can be shown that fronts 
with constant velocity appear only if the initial conditions for the particle densities decrease 
exponentially with the distance from the cathode. In this context, the compactly supported 
case is treated as an exponential decay with infinite argument, as was done in the previous 
subsection. 

The shielding factor formulation allows us to deduce the existence of accelerated fronts. 
This situation occurs when taking an initial ionization decaying at infinity slower than an 
exponential. For instance, by taking 



A 



as x — ► oo , 



with a positive. Near the front, we can take u close to 1, so that from equation 

du du A 
dt dx x a 



If we introduce an expression of the form 

u = FU 



x 



(51) 
we get 
(52) 

(53) 
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into equation we obtain 



ld_£ + tia -. )/aF d_F = _A (54) 



which can be approximated by 

^ dF A F (K\ 
a~di = T« ' (55) 

when a < 1 and for large t. Hence there might exist fronts whose position is located at the 
points 

' : C, (56) 



t l/a 

implying a superlinear propagation, i.e. an acceleration. This result illustrates the usefulness 
of our method, showing an unexpected behavior of the well studied planar fronts. 



IV. CURVED SYMMETRIES 



When the initial particle distributions or the initial electric field do not have planar 
symmetry, the ionization wave behaves quite differently from what we have described in the 
previous section. In particular, the amplitude and the velocity of the travelling wave are 
always not constant. The shielding factor formulation can be applied to those more general 
curved cases with few changes with respect to the planar case. We will use this formulation 
to treat the cases of cylindrical and spherical symmetries. In the cylindrical case, we will 
see that the velocity of the fronts varies in time as t -1 / 2 and the amplitude of the front goes 
as 1/t when the initial conditions for the charged particle densities decay sufficiently fast 
with the distance to the cathode. In the spherical case, the velocity goes as t~ 2 ^ 3 and the 
amplitude varies as 1/t as in the cylindrical case. Both cases can be dealt in a very similar 
way. 



A. Cylindrical symmetry 

First we analyze the case with cylindrical symmetry. We consider the experimental 
situation of two cylindrical plates with radius r and r\ ^> r , respectively. The space 
between the plates is filled, as in the planar case, with a non-attaching gas. A constant 
potential difference V is applied to the plates, so that V{ri) — V(r ) = V > 0. Then the 



14 



initial electric field E (r) between the plates is 

E ( r ) = --u r , B= J\ (57) 
r log(ri/r ) 

where B is a positive constant and r is the radial coordinate, ranging from r$ to r\ . An 
initial neutral seed of ionization po( r ) with cylindrical symmetry is taken, so that n e o(r) = 
n p o{r) = Po( r )- It is useful to change the spatial variable r to 

* = (58) 
so that equation (|34jl for the shielding factor w takes the form of the Burgers equation 



S+U5 = -»*(*)-«_/ « (59) 



Since the integral term in equation (|59|) depends explicitly on x, it is quite convenient to 
define the variable v(x,t) through 



r^/B/(2x) 

v(x,t) = / e~ 1/s ds. (60) 

Jy/B/(2x)u 

Now, as in the case of planar symmetry, we can integrate (|S9*jl along characteristics, trans- 
forming this equation into the system of ordinary differential equations 

I = <61) 

_ = -uv - po(x)u, (62) 

Given the definition of v in (|5Dji. by taking its time derivative, and using (JBTj) and (jS2j)) we 
close the above system with the equation 



fib = VB/2 -y/^/E/u -3/2 2 

dt 2 

^B]2 



'V^ x -^u (63) 



+ ^/Bj^e-V^^x- 1 / 2 [ uv + Po ( x )u] . 

Equations (JfiT|) . J02J), ffHTTjl . constitute a Lagrangian description of the problem. This dynam- 
ical system can be solved with appropriate initial conditions x(0) = xq, u(0) = 1, v(0) = 0, 
for any Xq, allowing us to obtain the profiles for the function u(x,t) at any time t. 
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FIG. 3: (Left) Electron density n e of an ionization wave with cylindrical symmetry for fixed time 
intervals vs coordinate x = r 2 /(2B), where r is the radial coordinate and B is a constant related 
to the initial electric field between the plates. A compactly supported neutral seed of ionization 
near the cathode (x = 0) has been considered as initial condition for the electron density. The 
solution converges into a shock front with decaying amplitude and constant velocity (in terms of 
the non-physical x coordinate). The dashed line is the analytical prediction for the amplitude of 
the shock front. (Right) Same conditions as in (left) but the electron density has been plotted vs 
the physical radial coordinate r. The velocity is clearly non-constant, as explained in the text. The 
dashed line is the analytical prediction for the amplitude. 




FIG. 4: Points are the positions r of the shock front shown in Fig. |3] vs time (in the inset, we plot 
the positions in terms of x = r 2 /(2B)). Solid lines are, in both cases, the analytical predictions for 
the position of the shock front as a function of time. 
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The solutions of the dynamical system given by equations (pjTj) . (JHH), (ESI), depend on the 
particular choice of the initial condition for the electron and positive ion densities. In these 
equations we have used a neutral initial condition given by n e0 (r) = n p0 (r) = po(r). Consider 
now the especial case in which the initial condition po( r ) f° r both densities is compactly 
supported, strictly vanishing beyond a certain point. One example of this behavior is given 
by a homogeneous thin layer of width 5 <C (r\ — ro) from r = r$ to r = ro + 5, i.e. 



By using this initial condition, we have plotted in Fig. El the electron density distribution n e 
corresponding to a given choice of the physical parameters Vo, Po, 8, r o and r\. The electron 
density has been calculated from the shielding factor u using the relation (|23|) and plotted as 
a function of x at fixed time intervals. In Fig. EKieft), the electron density has been plotted 
as a function of coordinate x. In Fig. fright), it has been plotted as a function of the 
radial coordinate r with the help of the relation ()58|) . What we can see from the figure is 
a shock with decaying amplitude, separating the region with charge and the region without 
charge. The numerical data allows us to measure the velocity of propagation of such front. 
In Fig. HJ we have plotted the position of the shock r/ as a function of time t. The velocity 
of propagation is clearly not constant. However, when we plot the position of the front in 
terms of x, one can observe the following linear relation (see inset Fig. 0]), 



which implies, in the original cylindrical variable r, an asymptotic behavior such that the 
position of the front depends on time as 




(64) 



Xf(t) =t + x , 



(65) 



r f (t) 



Vwi, 



(66) 



so that the velocity of the front behaves as 




(67) 



Using (JSTj) and (JBT?)) . this result can also be written as 



c r ~ E (r) 



(68) 
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showing a close similarity with the case (|4bj) of planar symmetry. Physically, it means that 
the shock front moves with the drift velocity of electrons as it should be expected. 

This behavior can also be obtained analytically, as well as the amplitude of the shock 
front. This is a considerable advantage of using the formulation in terms of the shielding 
factor. To do such computation, it is useful to write, locally near the front, the solution for 
u as 

u(x,t) = l-a(tMO, Z = x-x f (t). (69) 
This expression can be substituted into equation ([59*]). The computation is simplified if we 
note that the integral term in (J59)) is very small when i>1. We get 

a(*V(0 - a(ty(£)4(t) - a'(t)<p(£) + a 2 (*M£V (£) ~ °- ( 70 ) 
The only way this equation can be satisfied is by choosing 

x f (t) = t + x , (71) 

o(i) = (ifij- (T2) 
<pd ) = - 1 - x„), (73) 

where (3 is an arbitrary constant depending on initial conditions. Equation ()71j) is an ana- 
lytical proof of the numerical law (J(i5|) obtained for the position of the front in terms of time. 
To obtain the amplitude of the shock front, we use relation PHJ) to compute the electron 
density from the asymptotic solution What we get is 

{ 1 l+(x-t-x )/(t+t ) T<f , T 
t+t0 ^ {x ^ xo)/(t+to) ,x_z + x ^ (74) 
, X > t + Xq 

which implies that the amplitude of the front decays with time as 

n e (x f (t),t) = (75) 
£ t to 

The analytical curve (|T5j) has been plotted as a dashed line in Figs-Efleft) and fright). The 
agreement with numerical data is seen to be excellent, especially for large times. 



B. Spherical symmetry 

The physical case in which the initial electric field and the initial particle densities have 
spherical symmetry shows close similarities with the cylindrical symmetry case. The shield- 
ing factor formulation for spherical symmetry has been used in to analyze a typical 
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corona discharge. In this example, we have two spherical plates with radius ro and r\ ^> ro, 
in which a potential difference V(ri) — V(r ) = Vo > is applied. Note that, in this 
case, r is the spherical radial coordinate. The initial seed of ionization is neutral so that 
n eo( r ) — n po{ r ) — Po{ r )i an d the initial electric field E (r) between the plates is 

Eo(r) = ~Ur, C = V ^^. (76) 
r A T\ — r 



Changing the spatial variable r to 



y»3 



*=3CT (7?) 



the evolution for the screening factor takes the form of the Burgers' equation 

i c W 3 

du du . , fwi 1/s , ,_ n . 

i+ « s = -, ft W-»/ i(6)l/ "*, (78) 

where po(x) is the initial distribution of charge. This equation, as in the case of cylindrical 
symmetry (|59j) . can be integrated along characteristics. The results are very similar to 
that of cylindrical symmetry shown above For the case of sufficiently localized initial 
conditions, when the initial electron density strictly vanishes beyond a certain point, there 
appears a sharp shock with decaying amplitude, separating the region with charge and the 
region without charge. The velocity of propagation of such front is given by the relation 
between the position of the front and time: Xf(t) = t + xq. This implies, in terms of the 
original variable r, an asymptotic behavior 

r f {t) ~ (3C) 1/3 1 1/3 , (79) 

for the position of the front. The velocity of the front is then 

Cr (t) ^(3cy/ 3 r 2 / 3 , (80) 

or, in terms of the initial electric field (|76|). 

c, ~ Eo(r). (81) 

The analytical computation of the amplitude and propagation velocity of the shock can be 
done, by taking the shielding factor near the front as u(x,t) = 1 — a(t)</?(£), which gives 
exactly the same equation found in the cylindrical case (J7UJ). The reason of this is that 
the integral term is neglected in both cases. However, note that the relation between the 
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coordinate x and the physical radial coordinate is different in each case. The position of 
the front then satisfies Xf(t) = t + xq and the amplitude of the electron density front also 
decays with the law n e (xf(t),t) = l/(t + 1 ). Details and figures showing these results can 
be found in [15| . 

V. GEOMETRICAL DIFFUSION AND SELF-SIMILAR BEHAVIOR 

In previous sections we have found that, in the framework of the minimal streamer model 
without diffusion, when an initial seed of ionization is placed near the cathode, a travelling 
ionization wave develops towards the anode. The shape and the velocity of this wave depends 
on the asymptotic behavior of the initial particle (electron and positive ion) density. If the 
initial particle density is compactly supported, i.e. it vanishes beyond a certain point, then 
the travelling wave is a shock front, whose velocity is equal to the drift velocity of electrons, 
i.e. equal to the modulus Eq of the initial electric field. This behavior is found in the case in 
which the physical situation has planar, cylindrical or spherical symmetry. In the last two 
cases, the velocity, as the initial electric field, is not uniform. With respect to the amplitude 
of the electron density of the shock, in the planar case, it is constant during the evolution, 
but it decays as 1/t in the curved case (cylindrical or spherical symmetry). 

When the initial particle density is not compactly supported, but it does decay exponen- 
tially with the distance from the cathode, the shock front does not appear. In the planar 
case, we have seen that the velocity of the front is then constant and larger than the drift 
velocity of the electrons, and the amplitude is constant. However, if the initial particle 
density decay slower than an exponential with the distance from the cathode, accelerated 
fronts might appear. 

Now we are going to investigate the especial features that appear in a case with cylindrical 
or spherical symmetry when the initial particle densities are not compactly supported but 
they decay exponentially with the distance from the cathode. We will see that (i) a shock 
front does not appear, (ii) but a front with an asymptotic self-similar behavior, (iii) whose 
velocity depends on time in a similar way as the velocity of the shock front seen in the 
previous section. As a remarkable fact, we will note the appearance of a new type of 
diffusion effect, due to the geometry of the initial physical situation. 

Consider the case of an initial electric field with a cylindrical symmetry as (JH7J) . The 
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FIG. 5: (Left) Electron density n e of an ionization wave with cylindrical symmetry for fixed time 
intervals vs coordinate x = r 2 /(2B), where r is the radial coordinate. An initial neutral seed of 
ionization near the cathode (x = 0) has been considered such that it exponentially decays with the 
distance from the cathode. The solution converges into a spreading front with constant velocity in 
terms of x. (Right) Same conditions as in (left) but the electron density has been plotted vs the 
physical radial coordinate r. The velocity is clearly non-constant in terms of r. 

initial electron and positive ion densities are equal and have cylindrical symmetry, so that 
n eo( r ) — n po{ r ) — Po( r ) ; where r is the radial coordinate. We use the variable x = r 2 /(2B) 
as in the previous section. Now we take an initial neutral charge distribution that is not 
compactly localized, for example given by 

n e0 (x) = n p0 (x) = po(x) ~ e~ Xx , x > 1. (82) 

We can solve this problem numerically, integrating along characteristics the dynamical sys- 
tem given by equations (joTj) . (JrJ2*|) and (JcJH|) . The obtained electron density appears in 
Fig. shown in constant time intervals, proving that the shock front does not appear. In 
Fig. Efleft), we plot the electron density vs x for different times. What we see is a travelling 
wave with increasing thickness and, remarkably, the center of this front moves with constant 
velocity in the coordinate x in the same way as the shock front does. This means that the 
velocity of the front center has a similar behavior c r ~ t~ 1 / 2 that the velocity of the shock 
front that we analyzed in the previous section. In Fig. fright), we plot the electron density 
vs the radial coordinate r. 

By using the shielding factor formulation, we can prove that the asymptotic local behavior 
of the electron density near the front of the travelling wave is self-similar. In order to show 
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this property, we introduce 

u{x,t) = l-g{x,t), (83) 

in the equation ()59|) describing the evolution of the shielding factor. As we are analyzing 
the asymptotic behavior, we can take x ^> 1 near the front. Keeping the main order terms 
we obtain the equation 

| + d-,)| = 0, (84) 

which can be written as 

¥-^ = °' (85) 
if we define £ = x — t. Equation (}85|) is a Burgers equation whose solution g is such that it 
is constant along the curves given by 

I- < 86 » 

If g varies slowly in some region of size 5 X in space, one can consider g — G at that region, 
G being a constant. An intermediate asymptotic regime is then obtained, such that it is 
constant along £ = Gt and therefore we can assume 

g( x ,t)^g(^-). (87) 



Hence the asymptotic behavior of the electron density is given by 

Idu 1 £/(5 x t) d 



n,. 



9 (WW). (88) 



udt tl+g(£/(6 x t))d? 

Consequently, the asymptotic local behavior of the electron density near the front is self- 
similar, given by 

n e (x,t)~h(4i), (89) 



t \s x t, 

in which £ = x — t, and / is some universal self-similar profile. Hence the front presents a 
typical thickness given by 

£ c ^ S x t. (90) 

This result can be seen in Fig. |BJ in which the self-similar character of the asymptotic local 
behavior is shown. The consequence of this result is clear: even neglecting diffusion, the front 
spreads out linearly in time when the initial condition for the particle densities decreases 
exponentially with the distance from the cathode. This is a new and remarkable feature, 
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FIG. 6: Asymptotic behavior of the electron density shown in Fig. |3 We plot the quantity n e t vs 
£/(8\t), and a self-similar behavior is apparent. The front spreads in time, showing a new type of 
diffusion effect completely due to the geometry of the initial field distribution. 

first considered in and explained here, of the curved geometry. It is a diffusive behavior 
of the solutions of the minimal streamer model caused by the geometry of the electric field. 
It has been termed geometrical diffusion. 

In the spherical case, the same behavior is found, with the only difference being that 
the x coordinate is related to the radial coordinate r in a different way. We conclude than 
geometrical diffusion is a universal behavior. 



VI. CONCLUSIONS 



The conclusions of this paper are the following. We have made a thorough study of the 
properties and structure of anode-directed ionization fronts without diffusion, based on a 
minimal streamer model for non-attaching gases. This model includes impact ionization 
processes as source terms. The role played by the condition that the magnetic effects in the 
streamer discharges are neglected has been discussed. As a consequence, it has been proved 
that an electric shielding factor can be defined, and the physical quantities can be expressed 
as a function of it. 

A Burgers type equation is obtained for the evolution of the electric shielding factor. 
Thus, the analytical and numerical study of the ionization fronts can be performed by using 
a Lagrangian formulation. The power of this formulation makes it easier to treat the cases 
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of the non-homogeneous initial electric field in electric discharges with curved symmetries. 

We have applied this new formulation to the discharge between planar and curved elec- 
trodes (with cylindrical and spherical symmetries). When an initial seed of ionization is 
placed near the cathode, a travelling wave develops towards the anode. The shape and the 
velocity of this wave depends on the asymptotic behavior of the initial charged particle den- 
sities. If the initial density is compactly supported, then the travelling wave is a shock front, 
whose velocity is equal to the drift velocity of electrons. This behavior had been predicted 
for the planar case [6j], but we have found that a similar situation takes place in the case 
that the physical situation has cylindrical or spherical symmetries. We have derived power 
laws for the velocity and the amplitude of the shock fronts in the cases of curved symmetry. 
When we have cylindrical symmetry, the velocity of the shock front behaves as t -1 / 2 , and 
the amplitude behaves as t _1 . In the spherical case, the velocity of the shock front behaves 
as t~ 2 / 3 and the amplitude goes as t~ x . 

When the initial particle density is not compactly supported, but decays exponentially 
with the distance from the cathode, the shock front does not appear. In the planar case, we 
have seen that the velocity of the front is then constant and larger than the drift velocity of 
the electrons, as predicted in |6|. However, if the initial particle density decay slower than 
an exponential with the distance from the cathode, accelerated fronts appear. 

In the cases in which the physical situation has cylindrical or spherical symmetries and 
the initial ionization seed decays exponentially fast to infinity, we have seen that the velocity 
follows the same power laws as the compactly supported case. However, the structure of 
the travelling wave is rather different. We have proved that the asymptotic behavior of the 
charged particle densities is self-similar. Even if the diffusion has not been considered, the 
front spreads out linearly in time. This is a remarkable feature, first considered in |l5] and 
explained here. It is a diffusive behavior of the solutions of the minimal streamer model 
caused by the geometry of the electric field, and we have called it geometrical diffusion. 

Our analysis opens the way to consider geometrical effects in the stability of ionization 
fronts. 
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